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Abstract 



We use ISAJET to perform a detailed study of the missing transverse 
energy I$ T plus multi-jet signal expected from superparticle production at the 
CERN LHC. Our analysis is performed within the framework of the minimal 
supergravity model with gauge coupling unification and radiative electroweak 
symmetry breaking. We delineate the region of parameter space where the 
]$ T supersymmetry signal should be observable at the LHC and compare it to 
the regions explorable via searches for sleptons and for chargino/neutralino 
production. We confirm that, given a data sample of 10 fb _1 , nig ~ 1300 
GeV can be explored if nig ~S> nig, while nig ~ 2000 GeV can be probed 
if nig ~ nig. We further examine what information can be gleaned from 
scrutinizing this event sample. For instance, the multi-jet multiplicity yields 
information on whether squark production makes a significant contribution 
to the observed ~f$ T sample. Furthermore, reconstructing hemispheric masses 
may yield a measure of nig to ~ 15 — 25%. Finally, for favourable ranges of 
parameters, by reconstructing masses of tagged bb jet pairs, it may be possible 
to detect Higgs bosons produced via sparticle cascade decay chains. 
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I. INTRODUCTION 



A. Motivation and Framework 

The existence of weak-scale supersymmetry (SUSY), stabilizing the electroweak symme- 
try breaking sector of the Standard Model (SM), is a tantalizing hypothesis 0. This hy- 
pothesis has received some support from the observation that the running gauge couplings 
unify at a scale Mx ~ 10 16 GeV provided that the superpartners have masses ~ 1 TeV, 
while this unification does not occur in the SM 0. Hence, the search for supersymmetric 
particles is a high priority task for collider experiments. Expectations for super-particle 
masses are typically in the 100-1000 GeV range, whereas current limits on them are gener- 
ally below Mz/2 ||, although the CDF and DO experiments at the Tevatron probe gluinos 
and squarks as heavy as 150-200 GeV [HI. In the near future, LEP II will probe slepton 
and chargino masses up to mj, ~ 90 GeV, while Tevatron experiments, for favourable 
parameter ranges, may indirectly reach m g ~ 500 GeV in the main injector era, via the 
clean tri-lepton signal from W\Zi production These experiments have a good chance 

to discover weak scale supersymmetry, but they cannot exclude it. A thorough search for 
low-energy supersymmetry can only be made at the recently approved CERN Large Hadron 
Collider (LHC) or, alternatively, at e + e~ linear colliders with center of mass energies ex- 
ceeding 500-1000 GeV g^. 

Most early analyses were carried out within the framework of the Minimal Supersymmet- 
ric Model (MSSM), which is the supersymmetric extension of the SM with the least number 
of additional new particles and interactions necessary for phenomenology. As a consequence, 
i?-parity is exactly conserved, so that the lightest SUSY particle (LSP) is stable. Supersym- 
metry breaking is parametrized by the introduction of soft SUSY breaking mass terms and 
interactions consistent with the SM gauge group. Since each SU(3) x SU(2) x U{1) multi- 
plet of sfermions and gauginos has an independent mass, the most general parametrization 
of SUSY breaking requires a large number of mass parameters, and also many arbitrary 
trilinear scalar couplings. 

Without further assumptions about these parameters, any phenomenological analysis 
is essentially intractable. Motivated by grand unified theories (GUTs), we shall, following 
the early SUSY analyses in Ref. [PHTT| , further assume that the three DR gaugino masses 
originate in a common gaugino mass at some high scale, so that the ratios of these masses at 
the weak scale are proportional to the corresponding ratios of the squared gauge couplings. 
The 577(2) and U(l) gaugino masses are then determined in terms of m g . Motivated by 
supergravity models, these analyses further assumed a common physical mass for each of 
the light squarks, left and right sleptons, and sneutrinos. The Higgs sector of this model is 
determined at tree level by m# p , the mass of the neutral pseudoscalar, in addition to the 
parameters tan (3, the ratio of the two vacuum expectation values, and /i, the superpotential 
Higgsino mass (which also enter the gaugino- Higgsino sector). Finally, the weak-scale A- 
parameters, which mainly affect the phenomenology of third generation squarks were needed 
to completely specify the model. In what follows, and in the simulation program ISAJET 
jl2]| , we use the term MSSM to refer to this supergravity-inspired framework. 

In view of the fact that the additional assumptions for the MSSM are motivated by 
supergravity GUT models, it seems reasonable to seriously explore their implications. The 
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assumptions underlying these models, which differ in important ways from those of the 



MSSM, are outlined below. Several phenomenological analyses of supergravity models |13 
have recently appeared in the literature flI^]I^|I5|JI7| ■ In minimal supergravity with canon- 
ical kinetic energy terms, gravitational interactions communicate the effects of supersym- 
metry breaking in a hidden sector to the observable sector of particles. These gravitational 
interactions, being universal, lead to a common mass (mo) for all scalar fields and common 
trilinear and bilinear soft SUSY breaking (Aq and Bq) scalar interactions The resulting 
Lagrangian should be regarded as an effective Lagrangian of the theory below the Planck 
scale, with its parameters renormalized at some ultra-high energy scale close to Mpi anc ^. If 
we further assume that these interactions respect an (unbroken) GUT symmetry, the gaug- 
inos must all be part of a single multiplet, and so must have a common mass (mi/2) at the 
this scale. Similar considerations also apply to string models, which can lead to coupling 
constant unification without a grand unifying gauge group [jnj. The SUSY breaking sector 



of the model is then completely specified by these four parameters. 

For phenomenological analyses, the running parameters in the Lagrangian should be 
renormalized at the weak scale to sum the large logarithms arising from the disparity be- 
tween the high scale where the physics is simple, and the weak scale relevant to collider 
phenomenology. This is most conveniently done using renormalization group (RG) equa- 



tions [20] and taking the common values of the four high-scale input parameters as boundary 
conditions. The RG evolution splits the various masses and leads to a rich spectrum at the 
weak scale [fJTJ. The first two generations of squarks are approximately degenerate, so these 



models automatically satisfy constraints [£2| from the non-observation of flavour changing 
neutral currents in the kaon sector, while t\, the lighter t-squark and h\ ~ bi, the lighter 
6-squark, may be substantially lighter. Sleptons are lighter than squarks, and may be much 
lighter if gluinos and squarks of the first two generations have comparable masses. As a 
bonus, the RG evolution also gives the correct pattern of electroweak symmetry breaking 
2"B| , |2"0"[ for considerable ranges of input parameters. It is customary to eliminate Bq in favour 



of tan/? and to adjust the value of fi to reproduce the measured Z boson mass. Then the 
particle/sparticle masses and couplings are determined by the following parameters (along 
with m t ): 

m , mi/2, A Q , tan/5, and sgn(p). (1) 

In particular, the values of the MSSM parameters fi and mjj p are determined. Typically, \i 
is large, so that the lighter charginos and neutralinos are dominately gaugino-like, while the 
heavier ones are dominantly Higgsino-like. Also, m# p 3> My/ so the lightest Higgs scalar 
(Hi) resembles a SM Higgs boson with the other Higgs particles all relatively heavy. In what 
follows, we shall use the term SUGRA to mean supersymmetry with masses and mixings 
calculated from the parameters in equation (1). 



The simulation program ISAJET |L2] now allows the user either to specify the weak-scale 
MSSM parameters — the gluino, squark, stop, sbottom and slepton masses, A-parameters, 
fi, rriH p and tan (3 — or to use values of these parameters calculated from the SUGRA pa- 
rameter set (1). Because the SUGRA model is specified by just four new parameters, various 
experimental observables that can be determined in experiments at colliders become corre- 
lated and so provide non-trivial tests of the underlying assumptions ]T^JT7| ,[7[] . It is, however, 
worth remembering that SUGRA models are based on extrapolations of physics which may 



3 



ultimately prove incorrect. While these models are indeed very attractive, economical, and 
satisfy all phenomenological constraints, it may be worthwhile to test the sensitivity of model 
predictions by running ISAJET using the optional MSSM input set, which is a more gen- 
eral framework that can encompass models with non-universal soft-breaking sfermion mass 
terms. This may be particularly important for assessing the reach of future facilities. 



B. Phenomenological Overview 



Since strongly interacting sparticles are most copiously produced at hadron colliders, 
many of the early studies of supersymmetry at hadron supercolliders focussed on the detec- 
tion of gluinos and squarks P-4TT1. Recently, a number of papers addressing the detection of 
weakly interacting sparticles at the LHC have also appeared. In Ref. , it was shown that 
LHC experiments ought to be able to detect clean dileptons from slepton pair production if 
rrii ~ 200-250 GeV. Also, in Ref. [ 2~5\ , it was shown that LHC experiments ought to be able 
to detect the clean trilepton signal from W\Z 2 production over much of parameter space 
as long as the two body decays Z 2 — > H^Z\ and Z2 — > ZZ\ are kinematically inaccessible. 
Experiments at a high luminosity upgrade of the Tevatron may have a similar reach as the 
LHC if n is negative; for positive values of fi the branching fraction for the three body decays 
Z 2 — > UZ\ is strongly suppressed, and the corresponding reach at the Tevatron is somewhat 
smaller than at the LHC. 

Direct slepton pair and chargino/neutralino production takes place via weak interactions, 
whereas the strength of the LHC lies in the production of strongly interacting particles. In 
Fig. [I], we show the total production cross sections at -^i = 14 TeV for strongly interacting 
SUSY particles (gg, gq and qq), for charginos and neutralinos in association with squarks 
and gluinos, and for gaugino pairs {W\Z 2 and VP1W1), as a function of nig. We have assumed 
gaugino mass unification and 5 degenerate species of L- and R-squarks, taken /i = — nig and 
tan (3 = 2, and used the CTEQ2L parton distribution functions In Fig. |]a, we take rriq = 
trig, while in Fig. [l|& we take rriq = 2m g . In a), it is clear that strong sparticle production 
is the dominant production mechanism at the LHC for all values of m g from 300 GeV out 
to 2000 GeV. In b), strong sparticle pair production is dominant up to m g ~ 1100 GeV, 
after which chargino/neutralino production becomes dominant. In both cases, associated 
production mechanisms are sub-dominant cross sections, even though the combined mass of 
the produced sparticles can be much smaller than the mass of a pair of strongly interacting 
sparticles. The fact that strongly-produced sparticle pairs are the dominant production 
mechanism for squarks and gluinos as heavy as 1-2 TeV leads us to focus on signals from 
squarks and gluinos. 

Once gluinos and squarks are produced, they are expected to decay through various 
channels until the cascade terminates in the stable LSP (taken to be Z\) |27|]. These cas- 
cade decays lead to final states containing multiple jets, isolated and non-isolated leptons, 
and missing transverse energy (mainly from the undetected Z\§ and neutrinos). Rates for 
the various multi-lepton signatures have been presented in Ref. f28f; the most promising 



signatures appeared to be those containing same-sign isolated dileptons [p9|| , and isolated 



trileptons. The multi — jet + If T signal, while generally much larger than signals from 
multi-leptons, is less clean, due to irremovable backgrounds from various SM QCD induced 
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processes (multi-jet production, vector boson production in association with jets and heavy 
flavor production). 

Detailed studies of the missing energy plus multi-jet signal have been performed by the 
GEM JTTJI and SDC |K|] collaborations for the now defunct Superconducting Supercollider 



project |[30|| . More recently, detailed studies of this same signature have also been performed 
by the ATLAS Collaboration for the CERN LHC. In these LHC studies, performed 
within the framework of the supergravity-inspired MSSM, it was shown that values of m g ~ 
1300 GeV [m-g ~ 2000 GeV) could be probed for nig m g (m, ~ rn g ), given 10 fbT 1 of 
integrated luminosity. In addition, values of m g as low as 300 GeV could easily be probed at 
the LHC, so that there should be no "gap" in the explorable range of m g between Tevatron 
searches and future LHC searches. 

In the present paper, we further scrutinize the missing energy plus multi-jet signal. Our 
goals are multiple: 



1. We update previous results on multijet +If T cross sections by using the latest 
ISAJET 7.13 simulation for sparticle production and cascade decays. 

2. We evaluate the reach of the LHC via some set of optimized cuts, valid across a large 
range of sparticle mass choices. 

3. We present our reach results in the parameter space of the minimal supergravity model 
with gauge coupling unification and radiative electroweak symmetry breaking. The 
rather small parameter set yields a correlated spectrum of SUSY particle masses, and 
allows one to compare the regions of parameter space that different search experiments 
can probe. 

4. We examine what information can be gleaned from a sample of signal events in the 
missing energy plus multi-jet channel. For instance, can one tell whether the signal is 
mainly due to gluino production, or whether qq and qg production is also significant? 
Can one gain some sort of mass measurement for the gluino or squark? Can one 
identify the presence of other particles or sparticles by sifting through the debris of 
the cascade decay? 

The first three of these points are addressed in Sec. II of this paper, where we mainly map 
out regions of SUGRA parameter space explorable by LHC experiments using the multi-jet 
plus missing energy signature. We also compare these regions with those where various 
slepton and chargino/neutralino searches are expected to result in observable signals. In 
Sec. Ill, we address the questions raised above in item (4), and find that to some degree, 
and at least in some cases, all the questions raised can be answered in the affirmative. We 
summarize our results and present our conclusions in Sec. IV. 



II. REACH OF THE LHC IN SUGRA PARAMETER SPACE 

A. Event simulation 

In this paper, we work within the framework of the minimal SUGRA model, as imple- 



mented in the ISAJET 7.13 subprogram ISASUGRA [IT]. The SUGRA parameter set (1 
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can be directly entered into ISAJET, and ISASUGRA then calculates all the SUSY particle 
masses and mixings. Briefly, starting with the precisely measured gauge couplings at scale 
M z , ISASUGRA evolves the three MSSM gauge couplings and third generation Yukawa 
couplings to a mass scale (determined to be ~ 2 x 10 16 GeV) where cti and a 2 unify. At this 
approximate unification scale, a 3 is set equal to a\ and a 2 , and the GUT scale boundary 
condition values for mo, mi/2, A and B are implemented. From the approximate unifica- 
tion scale, the various soft SUSY breaking masses, gauge and Yukawa couplings are evolved 
via 26 RGE's to their weak scale values. We use 2-loop RGE's for gauge couplings, but only 
1-loop RGE's for soft breaking terms. Weak scale sparticle threshold effects are included 
within the gauge coupling evolution. Soft breaking masses are evolved only down to the 
scale value equal to their mass value. At the electroweak scale, the 1-loop effective potential 
is minimized, allowing the replacement of B by tan/3, and evaluating the magnitude (but 
not the sign) of ji in terms of M Z - This procedure is iterated until a stable solution for all 
sparticle masses is obtained, consistent with grand unification and electroweak symmetry 
breaking. 

After calculation of the spectra and couplings as detailed above, ISASUGRA then cal- 
culates all available sparticle decay modes and branching fractions. Currently, the program 
is valid only for tan/3 ^ 10; for very large values of tan/3, the effects from sbottom and 
stau mixing, not yet included, become important. It also assumes that the LSP is Z Xl which 
is not true for all choices of the parameters in set (1). Next, ISAJET generates all lowest 
order 2 — > 2 subprocesses for sparticle pair production according to their cross sections. The 
produced sparticles are then allowed to decay via the various possible cascades. Initial and 
final state parton showers are implemented, as well as hadronization and beam fragment 
evolution. 

For detector simulation at the LHC, we use the toy calorimeter simulation package ISA- 
PLT. We simulate calorimetry covering —5 < rj < 5 with cell size Arj x A</> = 0.05 x 0.05. 
We take the hadronic energy resolution to be 50%/V^ © 0.03 for \r]\ < 3, where © denotes 
addition in quadrature, and to be 100%/\^E © 0.07 for 3 < \r]\ < 5, to model the effective 
Pt resolution of the forward calorimeter including the effects of shower spreading. We take 
electromagnetic resolution to be 10%/ \/E © 0.01. Jets are found using fixed cones of size 
R = VAr/ 2 + A0 2 = 0.7 using the ISAJET routine GETJET. Clusters with E T > 100 GeV 
and |r/(jet)| < 3 are labeled as jets. Muons and electrons are classified as isolated if they 
have pt > 20 GeV, \f](£)\ < 2.5, and the visible activity within a cone of R = 0.3 about the 
lepton direction is less than .&r(cone) = 5 GeV. 

B. Signal versus background 

To evaluate signals from supersymmetry in the multi-jet +1$ T channel, we generate 
all possible supersymmetric subprocesses using ISAJET. Major physics backgrounds come 
from various SM processes which can give large amounts of missing transverse energy due 
to neutrinos produced in events, due to mis-measurement by calorimeter cells, and due to 
dead regions of the detector. We evaluate to following SM background processes: 

• QCD multi-jet production (e.g. gg — > gg etc., where extra jet activity comes from 
parton showers, (this includes heavy flavor bb and cc production), 
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• Z+ multi-jet production, where Z — > vv or rf, 

• W+ multi-jet production, where W — > £^ or rz/ T , 

• tt production and decay. 

We impose a series of cuts to extract signal from the vastly larger SM production cross 
sections. Since we wish to detect gluinos or squarks over a large mass range ~ 300-2000 GeV, 
different cuts are needed to optimize signal/background depending on the sparticle mass. 
After exploring a large variety of possible cuts for various SUGRA parameter choices, we 
arrived at the following requirements: 

• no isolated leptons, 

• transverse sphericity St > 0.2 to reduce QCD dijet background, 

• number of jets n(jets) > 2 (jets as defined above), 

• transverse plane angle A0(^ T , j c ) between 1$ T and closest jet is 30° < A(f> < 90°. 

After these mass-independent cuts, we apply a variable cut with the parameter E T chosen 
depending on the gluino and squark masses: 

• E T > E T and E T ( 3l ), E T { ]2 ) > E T . 

The results with these cuts for signal and total background levels are shown in Fig. |2|, 
versus the cut parameter E T . We show the signal cross sections for the six cases listed 
in Table I, which roughly correspond to nig ~ 300, 800 and 1300 GeV, where rrig ~ rrig 
or rriq = (1.5 — 1.7)rrig. We also show the estimated background from the sum of all SM 
processes. Since the total SM cross section corresponds to about 10 15 events, it is obviously 
neither technically possible nor physically reasonable to generate realistic statistics for it. 
Instead, the various SM processes have been generated in several overlapping p? ranges to 
obtain a reasonable estimate for all E T . The background is discussed in more detail in 
Sec. lie below. 

From Fig. [| we see that for E T = 100 GeV, case 1 and 2 with rrig ~ 300 GeV are 
easily visible at levels of 5-10 above background, while the cases with heavier gluino masses 
are well below background. As the cut E T is increased, the SM background decreases 
much more quickly than the heavy gluino signals. For E T > 150-200 GeV, the cases with 
rrig ~ 800 GeV begin exceeding background. To see the signal from rrig ~ 1300 GeV above 
background, a hard cut of E T > 300-400 GeV is needed. This agrees qualitatively with the 
results of Ref. |JT), where two sets of cuts (soft and hard) were advocated to see light or 



heavy gluinos. As we will see shortly, gluino and squark mass values across the range of 
rrig ~ 300 up to rrig ~ 1300-2000 GeV ought to be detectable. In particular, there should 
be no "undetectability window" in m g between Fermilab Tevatron experiments of the Main 
Injector era, and LHC experiments 

Fig. || also suggests a way to get a crude estimate of the gluino or squark mass. One 
plots the observed cross section versus E T (for gradual stiffening of cuts), and compares with 
background expectations. The approximate range where rrig might lie can be obtained by 
measuring event rates for several values of E T , starting at the point where observation begins 
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exceeding expectations from SM processes, and comparing against Monte Carlo expectation 
of the signal. 

In Fig. |2], we have fixed A = and tan/3 = 2. Variations in A Q mainly affect third 
generation squark and slepton masses, so our results are relatively insensitive to different A 
values as long as the light stop t\ is not driven to too small a value. As an example, we plot 
in Fig. ||a the cross section after cuts (taking = 300 GeV) for m = m\/2 = 300 GeV 
(case 3) and also for mo = 2m\/2 = 600 GeV, for five choices of Aq. Little variation is seen 
when A varies over the range shown. Multi-jet+^ T signals are also relatively insensitive to 
variations in tan f3. Case 3 is also plotted in Fig. for five choices of tan j3 and again shows 
only small differences in signal cross section. 

C. Reach of the LHC in SUGRA parameter space 

Our next goal is to evaluate the reach of the LHC via the multi — jet + If T signature in 
SUGRA parameter space. From Fig. |2j, it is clear that for very heavy gluinos, a large value 
of Ej, is desirable to enhance the signal relative to the background. For the heaviest gluinos 
and squarks observable at the LHC, which are roughly the heaviest consistent with SUSY 
relevant to electroweak symmetry breaking, we want to make a cut at about Ej> = 500 GeV. 
Only a very small fraction of the SM cross section passes such a cut, so it is necessary to 
combine events generated in several ranges of hard scattering pt (Pt S ) f° r each SM process 
that can produce backgrounds. Table II lists the ranges generated for each process, the 
total number of events for each range, and the number of Monte Carlo events and the 
corresponding cross section passing the cut E^ = 500 GeV. Fig. || shows the background 
cross sections for each process and range of p^ s vs. E^ for those values of E^ for which 
it is non-zero for the available statistics. While in retrospect a more uniform distribution 
of Monte Carlo events in log pt would have given a better estimate of the background, the 
samples listed in Table II are sufficient to provide an estimate of the SM background for all 
relevant values of E^. In particular, for each E^ there are enough events in some range for 
each process to give an estimated background cross section greater than that from lower p T 
ranges that produce no events. Therefore, we estimate the SM background as a function of 
E^ using the 95% upper limit calculated from the ranges giving a nonzero number of events, 
setting the background from the ranges giving zero events to zero. In particular, the total 
background from Table II then comes out to 1.86 ± 0.36 fb, which gives a 95% CL upper 
bound on the total BG to be a = 2.44 fb. We conservatively use this value to obtain our 5<r 
reach. 

Since ISAJET generates higher order QCD processes using the branching approximation, 
it does not necessarily give the correct background for SUSY signatures, which typically 
involve "round" events with jets far from the collinear limit in which the branching approx- 
imation is correct. Fortunately, we find that our estimates of the SM background is smaller 
than the SUSY signal, so the SUSY mass reach is not very sensitive to the background 
estimate. Also, the significance of any signal is not dependent on the calculation of the SM 
backgrounds, since these can be determined from other data. Backgrounds from neutrinos 
from b and c decays can be checked against the pt distribution of non-isolated muons. Back- 
grounds from Z — > vv can be checked against measurements of Z — > £ + £~. Backgrounds 
from W and ti can be checked against distributions of isolated single and double leptons. 
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Detector induced backgrounds can be checked using jet-jet and 7-jet events. Given this 
data, it should be possible to determine the SM background, including the multijet and St 
cuts, accurately and so assign a significance limited only by statistics to any excess from 
SUSY. 

The 5cr background level calculated as described above is shown by the horizontal line 
in Fig. [|. (Notice that the signal to background ratio is very close to unity even for a 
signal just at this level so that we do not need to impose any other S/B requirement for 
the observability of the signal.) On the same plot, we show signal rates versus m g for two 
cases: i) squares for m = m 1 / 2 (rrig ~ rrig), and ii) triangles for m = 4mu2 ( m q m g)- 
From this plot, we see that for rrig ~ m g , the LHC should be able to probe to m g ~ 2000 
GeV, while for rrig > rrig, LHC should be able to probe to nearly m g ~ 1300 GeV. These 



results are in remarkably close agreement with the reach calculated in Ref. |JT|, which used 
somewhat different MSSM parameter choices and different cuts. 

In Ref. [|I7]], it was shown that in the minimal SUGRA model, the m vs. m.1/2 plane 
forms a convenient panorama in which to plot various constraints. In Figs, f^^, we plot the 
reach of the LHC in the multi-jet+^ T channel, using the 5a constraint for 10 fb _1 of collider 
data. For convenience, in frame a) we show contours in SUGRA parameter space; in frame 
b) we show corresponding gluino and squark (averaged over the 4 first generation squarks) 
mass contours. All the figures take A = and m t = 170 GeV. Fig. |6| plots contours for 
tan (3 = 2 with fi < 0, while Fig. [7| plots for the same value of tan /3, but with fi > 0. Finally, 
Fig. H shows results for tan/3 = 10 and \i < 0. 

In these figures, the region shaded with bricks is excluded on theoretical grounds: either 
the correct electroweak symmetry breaking pattern is not obtained (or gives the wrong value 
of Mz), or the LSP is not the Z\, but is instead the W\, the f\, £r or v L . In addition, for some 
values of A , the light stop t\ can be driven tachyonic, so that electromagnetic and colour 
symmetries are spontaneously broken. The hatched regions are excluded by experiment. 
The experimentally excluded regions are formed [T7J by combining bounds from LEP |3| 
that rriyy > 47 GeV, > 60 GeV, and rriy > 43 GeV. In addition, the latest bound from 
CDF/DO on multi-jets+i^ events 0] is included. 



In Ref. |p4|| , it was shown that LHC ought to be able to probe sleptons with mass 
mj ~ 200-250 GeV. The area to the left of the contour labelled £(200) denotes this region in 
SUGRA parameter space. Also, it was shown in Ref. that (for 11 < 0) LHC experiments 
ought to be able to explore much of the region below the contours labelled Z2 — > Z\Hi and 
Z 2 — > Z\Z (above which the so-called spoiler modes become accessible), by searching for 
the clean trilepton signature of W\Z 2 — > 3£. This region also corresponds approximately 
to the maximal reach of Fermilab Tevatron collider experiments (although large holes of 
non-observability exist for some ranges of parameters due to a strong suppression of the 
Z 2 — > Z\tl branching fraction, especially for positive values of 11). Furthermore, LEP II 
should be able to explore the regions below the contours marked if; (90) and Wi (90) where 
the lightest Higgs boson and the chargino are, respectively, lighter than 90 GeV. 

In each of Figs. ^]a-|^a, we see the 5a 1$ T reach is maximal for m ~ | m i/2 ~ 500 
GeV, and drops to intercept the theoretically excluded (due to m fl < m^) region around 
m ~ 200 GeV; this fall-off is due to the fact that as m decreases, sleptons become very 
light, resulting in the presence of many isolated leptons in the final state of gluino and squark 
cascade decays. The isolated lepton veto thus suppresses somewhat the multi-jet cross 
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section as rriQ decreases; we may, however, expect that the multilepton cross section is large 
in this range of m . For not too large values of m , squark masses are somewhat lighter 
than (or comparable in mass to) gluinos, so that the combined cross section for gg, gq and 
qq production is very large. In this region, the maximal reach of the LHC is obtained: 
rrig ~ 2000 GeV can be explored. For mo very large compared to mi/2, squarks and sleptons 
are far heavier than gluinos. In this case, superparticle pair production is dominated by 
just gg production (for m g < 1100 GeV). Nevertheless, a mass reach of m g ~ 1300 GeV is 
obtained. 

Recently, upper limits on sparticle masses have been obtained by requiring no fine-tuning 
in minimal SUGRA models J32|. These limits, which are somewhat subjective, suggest 
m>g, rriq ~ 700-800 GeV. Comparison of these values with the calculated reach of the LHC 
suggests that LHC can perform a complete scan over the parameter space of minimal SUGRA 
models. 



III. CHARACTERISTICS OF THE MULTI-JET +# T SIGNAL 

If supersymmetry is discovered, it will be important to see what information can be 
gleaned from signal events about the properties of super-particles. We have already men- 
tioned that a rough determination of m g and rrig can be made by noting the size of the 
signal cross section above background for different choices of the cut parameter E^. In this 
section, we focus on several examples of intrinsic event properties, and the information they 
provide. 



A. Jet multiplicity: Detecting squarks in the Jf, T sample 

If m g is sufficiently smaller than rrig, then gluino pair production dominates, and each 
gluino decays typically via a lengthy cascade resulting in numerous final state jets, e.g., 
g — > qqZi, g — > qqWi or g — > tt\, followed by further Z i; W iy t and t\ decays. In contrast, 
if rrig < rrig, squarks decay through a more abbreviated cascade, via qi — > qWi, qZ iy and 
q~n — > qZi (Third generation squarks can, however, decay to charginos via their couplings 
to the Higgsino component). For right squarks, qn — > qZ\ dominates over large regions 
of parameter space. Since hard jets are most likely to come from the primary decay, in 
multi-jet +^/ T events, for a fixed value of rrig, one frequently expects a higher jet multiplicity 
if rrig > rrig, and gg production dominates, than if rrig < m g , and qq and qg production 
dominates. (An alternative method involving the study of the charge asymmetry in the 
same sign dilepton sample was suggested in Ref. p8| .) 

This is explicitly illustrated in Fig. |^ where we plot the jet multiplicity after cuts (for 
the Ej, value listed) for each of the SUGRA cases listed in Table 1. We have relaxed the 
jet cut to require only pr(jet) > 50 GeV for this plot, so that we increase the sensitivity to 
lower energy jets produced further down the cascade decay chain. The histograms shown 
include both signal and background contributions. In Fig. |9]a, we show jet multiplicity for 
the two rrig ~ 300 GeV cases (with E? = 150 GeV), where case 1 has rrig ~ rrig, and case 
2 has rrig ^> m g . We see that case 2 is dominated by 4 and 5 jet production, indicative of 
gg production and decay, whereas case 1 has a significantly lower average jet multiplicity, 
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and produces mainly 3 and 4 jet events more characteristic of a large component of qq 
production. In Fig. the two cases for m g ~ 800 GeV again show a larger jet multiplicity 
for the case where m g <C m§. Finally, Fig. |9]c shows the jet multiplicity for the two cases 
with nig ~ 1400 GeV, and again the m g <C rrig case has larger jet multiplicity. For case 6, 
we see that the jet multiplicity is large for n(jet) = 3, and then falls for n(jet) = 4, and rises 
again to a maximum for n(jet) = 5. This is due to a large background contribution in the 
n(jet) = 3 bin. 

In all cases, the n(jet) distribution has a wide range of smearing. This is due to multiple 
hard partons produced at various stages along the cascade decay chain, but is also due to 
substantial initial and final state QCD radiation. Nonetheless, the final distributions do 
show that, if we have some idea of the gluino mass, the measured jet multiplicity can give a 
handle on whether or not the event sample contains a significant q production component, 
thus reflecting the relative values of m g and rrig. Although we have illustrated this only for a 
few cases, we expect this to be a generic feature. The multiplicity in gluino events can only 
be reduced if the branching fraction for the radiative decay g — > gZ% becomes very large. 

We have also examined several other techniques for detecting the presence of squarks in 
the Ifirp sample. These are based on the expectation that heavy squark events (which form 
only a small fraction of the total sample) are expected to be more spherical and have larger 
jet multiplicities. For the first two cases in Table I, we have examined (i) the transverse 
sphericity distribution, (ii) the scatter plot of the transverse sphericity (St) versus the 
"bigness" B = |^ T | + S|pr(jet)|, (iii) scatter plot of B versus ^(max), the largest transverse 
momentum relative to the sphericity axis and (iv) the scatter plot of B versus n(jet). Only 
the last of these distributions appear to show any significant difference. Because this last 
distribution is correlated with the jet multiplicity distributions that we have already studied, 
we do not show it here. 



B. Gluino mass measurement 

Although measurement of sparticle masses can proceed with significant precision at e + e~ 
colliders |T3| (given sufficient energy and luminosity), the situation is expected to be much 
more difficult in the environment of a hadron supercollider. This is especially true of mea- 
suring the mass of the gluino at the LHC Even in the ideal case where all final state decay 
products of a gluino are tagged and precisely measured, the invariant mass distribution 
created will be a broad distribution between zero and m g — . Real events will contain 
overlapping decay debris from both of the sparticles produced in the subprocess, in addition 
to significant amounts of QCD radiation, and imperfect detector effects. 

In the past, various methods have been proposed for sparticle mass measurements at the 
LHC: 



In Ref. [Q, associated production processes such as gZ\ were examined using parton 
level event generators. If these events could be singled out, then the ambiguities from 
producing and decaying two strongly interacting sparticles are by-passed. In practice, 
rather hard cuts were required to separate the gZi events from gg events. Ultimately, 
it was concluded that this reaction might be of use only if m g < 350 GeV. 
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In Ref. , it was shown that — m-^ could be measured with significant precision 

from the end-point of the dilepton invariant mass distribution. However, the W\Z 2 — > 
3£ signal upon which this is based is only observable in a limited region of parameter 
space. 



• In Ref. p9fl , it was claimed that m g could be measured to 15% by focussing on same-sign 
isolated dilepton events from gg production, where each gluino decays via g — > qqWi, 
with W\ — > lvgZ\. However, these calculations may be overly optimistic, since this 
study considered only a single production mechanism (gg production- whereas same 
sign dileptons can come from various SUSY sources), a single cascade decay chain, and 
neglected effects of additional QCD radiation. 

Here, we seek to measure the gluino mass in multi-jet+i^p events, following a similar 
path to Ref. [p9[| . We proceed as follows. First, we have a rough estimate of m g and 



rrtq by examining signal to background levels versus cut parameter E T . After obtaining a 
relatively clean sample of signal events for an appropriate E T value, we divide the event 
into two halves, using eigenvectors from constructing the transverse sphericity St- We then 
construct the invariant mass of jets in each of the two hemispheres, using only jets with 
Pt > E T . Events are rejected if there is only one qualifying jet in each hemisphere. We 
then plot as M est the maximal value of the mass of each hemisphere. Results are shown in 



Fig. [10] using m = m 1 / 2 , A = 0, tan/3 = 2 and /i < 0. The corresponding value of m g is 
shown in each frame. Background level is shown by the shaded regions. In these plots, as 
usual, all supersymmetric subprocesses are contributing, with full cascade decays and QCD 
radiation effects, and detector smearing as given in Sec. Ila. 

We see in Fig. |T0|a, using E T = 100 GeV for m g ~ 300 GeV, that the M est distribution 



does indeed have a large smear due to the various effects listed previously. However, at 
least with our toy detector simulation, the two values of m g , which differ by ~ 15%, appear 



distinguishable. We plot values of M est for m g ~ 800 GeV and E T = 350 GeV in Fig. |10|fc, 



and plot for m g ~ 1500 GeV with E T = 600 GeV in Fig. [TTJc. Again, values of m g differing 



by ~ 15% appear distinguishable. We note, however, that using a too small value of E T , or 
too small a value of pr(jet), can lead to a large amounts of smearing in these distributions, 
and loss of distinguishability. These distributions appear workable only for rather narrow 
ranges of cuts that guarantee sufficient hemispheric separation of event debris, and that only 
the hardest debris is used in constructing M est , so that contamination of M est is avoided. 
In Fig. [TTJ, we show the same M est plot, but now for m = 4mi/ 2 , so that m g ^> m g 
(we show only the first two cases in Fig. |Tt] since the last case is not observable above 
background). The distributions appear somewhat less distinguishable than in Fig. |1(], so 
that a mass resolution of 15% may be more difficult to attain in this case; nevertheless, 
distributions from gluino masses differing by 25% appear to be distinguishable. The jet 
multiplicity and/or 5-multiplicity (see next section) may have to be used to distinguish 
whether trig ^> m g or m g ~ m g . 

C. Detecting Higgs bosons via i?-jets in cascade decay debris 

Multi-jet +1$ T events from gluino and squark cascade decays are expected to be unusu- 
ally rich in heavy flavor content (mainly B mesons, and t-quarks, if kinematically allowed) 
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compared to SM processes [|9|, |15| , |3~6jj . This is due to a number of effects. 

• The ii and hi masses are driven to values lower than the other squarks. This means 
decays like g — > tti or g — > bbi can be kinematically allowed, and dominate decay rates 
for large enough values of m~ g . Even if such decay modes are not allowed, three-body 
g decays to tiZi, tbWi and bbZi are enhanced relative to decays to other squarks due 
to the larger propagator factor for the lighter top and bottom squarks. 

• The various two- and three-body g decays to third generation quarks and squarks are 
further enhanced by the large third generation Yukawa couplings ||34[ . 



• Zi — > Zjbb decays are also enhanced by Higgs boson mediated decays |T3| . 

• Higgs bosons can be produced in large quantities in g and q cascade decays. The 
various Higgs bosons frequently have dominant decays into 3rd generation particles, 
enhancing the heavy flavor content of g and q events |J. 

In Table III, we list the average 5-meson multiplicity in multi-jet+.S T events, for nine 
cases of SUGRA parameters. To construct our tabulation, we have required 1$ T > 200 
GeV, S T > 0.2, and at least two jets with pr(jet) > 100 GeV and |r?(jet)| < 2.5. We then 
examine all jets with > 50 GeV in the same rapidity interval. If there is a 5-hadron 
within AR = 0.5 of the jet, it is tagged as a B with 50% probability; otherwise, it is (mis)- 
tagged as a B with 2% probability |3"T |. From Table III, we see that in cases 1 and 2, for 



mo = 100 GeV, the B multiplicity is rather low, < % >~ 0.9, compared with cases 3-9, 
with higher values of m . This is because all squarks and sleptons are relatively light, so 
any enhancements in production and/or decay to third generation quarks are small. For 
a fixed value of rrig, the highest B multiplicity occurs for the cases with large mo, which 
are dominated by gg production, with enhanced cascade decays to heavy flavors: here, 

< ub >~ 1.6. Cases 3-5, with intermediate vaules of mo, also have intermediate values of 

< ub >~ 1.1. We also see that, as discussed in Sec. Ilia, the average jet multiplicity is 
smaller for m small, and larger for mo large. 

Can one see direct evidence for Higgs production in cascade decay events? The light 
Higgs Hi is usually produced somewhere down the cascade decay chain via Z 2 — > Z\Hi. It 
has been shown in Ref. |35| that only in very limited regions of parameter space is the decay 



Hi — > 77 likely to be visible. We have investigated cases 1-9 of Table II to see whether the 
Hi — > bb decay can also be observed in cascade decay events. We do this by constructing the 
invariant mass m h i of all tagged -B-jet pairs. Results are shown in Fig. First we show 



case 1 in Fig. |12|a, where m = 100 GeV, and Z 2 dominantly decays to real sleptons and 
sneutrinos, so that Higgs production in SUSY events is minimal. In this case, a Higgs peak 
is not necessarily expected; we see a continuum distribution with only a slight bump at rriH e - 
In Fig. [T2]& (case 3), where both gluino and squark pair production is important, Z 2 — > Z\Hi 
occurs at 94% branching ratio, the Higgs bump at m b i = 89 GeV stands out against the 
continuum from other tagged -B-jet pairs in the SUSY events. This is despite the fact that 
decays to b and t squarks constitute 60% of the gluino decays. Higgs production from 
the cascade decays of the first two generations of squarks remains observable. Finally, in 



Fig. pic, we show the distribution for case 8, where Z 2 — > Z\Hn occurs with 98% branching 



ratio. However, in this case there are very many other sources of tagged .B-jet s from g 
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cascade decays, especially g — > tiZx^, bbZ\^ and g — > tbWi. Moreover, since gg production 
is by far the dominant process, there is a very high B multiplicity; then, the Higgs boson 
is frequently produced along with other £?-jets, so that the combinatorial background tends 
to wash out the Higgs bump. The Higgs signal was also detectable at the appropriate Higgs 
mass in case 4; the signal did not stand out for the other four cases in Table III. 

In conclusion, it appears that while the 5-multiplicity is frequently enhanced above SM 
expectations, the Hi — > bb signal can be identified only in limited regions of the parameter 
space: obviously, we need a substantial number of tagged B-jets from Higgs decay, but also 
that these events should not (as, for instance, in case 3) simultaneously contain other tagged 
5-jets which would then produce a large combinatorial background. Another requirement 
for the detection of the Higgs boson bump is that the rate from other sources of tagged 
5-jets (e.g. g — > tit and bbi decays, where the squarks may be real or virtual) should not 
be overwhelmingly large. 



IV. CONCLUSIONS 

In this paper, we have performed a detailed analysis of the multi-jet+^ T signal expected 
from production and decay of supersymmetric particles within the framework of the mini- 
mal supergravity model, in which sparticle masses and couplings are fixed in terms of the 
parameter set (1). Assuming an integrated luminosity of 10 fb - , we have shown that the 
Ifirp signal (mainly from gluino and squark production and decay) should be observable in ex- 
periments at the LHC for gluino masses ranging from values accessible to Fermilab Tevatron 
experiments, up to nearly m g ~ 1300 GeV (m g ~ 2000 GeV ) for rriq 3> m g (for rriq ~ Trig), 
confirming earlier studies by the ATLAS collaboration |3~I|. We expect that the reach (in 



terms of sparticle masses) is not very sensitive to the details of the model as long as R-parity 
is conserved. Comparing the LHC reach in terms of m g and rriq with somewhat subjective 
upper bounds from fine tuning constraints (recall that the resolution of the fine-tuning 



problem provided the original motivation for phenomenological SUSY), it seems likely that 
LHC can perform a thorough search for i?-conserving weak scale supersymmetry. 

We show the SUSY reach via the multi-jets+^ T channel in the m vs. plane of the 
minimal SUGRA model with gauge coupling unification and radiative electroweak symmetry 
breaking in Fig. The relatively small parameter set yields a complete, highly correlated 
sparticle spectrum, so that the plot of the reach in the I$ T channel can be compared with 
previous studies on the reach for sleptons and charginos/neutralinos. By comparing these 
various regions in Figs. |]-|||, we see that of the channels studied, the reach in multi-jets+-S T 
is by far the largest. However, it should be kept in mind that multilepton signals (e.g. 
same-sign dileptons and trileptons) from cascade decays of gluinos and squarks may also 
probe some or all of this region. These signals would be especially important if mo is very 
small, or if gluino decays to third generations fermions and sfermions are strongly enhanced. 
These multilepton signals could then provide a striking confirmation of a supersymmetric 
signal discovery in the multi-jet+^ T channel. 

We have also studied what further information about sparticle properties can be obtained 
by a careful study of the If T sample. We have used ISAJET to demonstrate that jet 
multiplicity can be a useful tool for detecting the presence of squarks in the tf T sample, 
and indirectly inferring the squark to gluino mass ratio. Of course, the distribution of jet 
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multiplicity is sensitive to the sparticle mass, so that this method is useful only after an 
estimate of the mass is obtained. We have further demonstrated that it might be possible 
to obtain a measure of the gluino mass by dividing each event in two halves, and using 
the greater of the masses in the two hemispheres as an estimator of rrig. We see from the 



distributions in Figs. |10| and |TT] that gluino masses differing by 15-25% might be possible 
to distinguish. Since we have included the production of all sparticles in our simulation 
as well as contamination to the hemispheres from QCD jets, we believe that it would be 
worth testing this technique to see if it survives a detailed detector simulation. Finally, 
we have studied if signals from Higgs bosons produced by cascade decays of gluinos and 
squarks and decaying via H e — > bb are detectable. We conclude that, with reasonable B 
tagging efficiencies, this is possible but only for favourable ranges of parameters, where 
events containing Higgs bosons are relatively free of other B-jets, and further, that S-jets 
from other SUSY sources do not overwhelm the Higgs signal. It should nonetheless be kept 
in mind that SUSY events frequently tend to be rich in central B-jets, so that this may 
provide another handle for distinguishing SUSY from the SM, and gaining information on 
the squark to gluino mass ratio. 
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TABLES 



TABLE I. Six parameter choices used for cross section evaluation of multi-jet +I$ T events. We 
take A = 0, tan/3 = 2 and sgn(/i) < 0. Also, m t = 170 GeV. 



Case tuq mi/2 m g m q 

1 100 100 290 270 

2 400 100 310 460 

3 300 300 770 720 

4 1200 300 830 1350 

5 600 600 1400 1300 

6 2000 500 1300 2200 



TABLE II. Details of background calculation after cuts using cut parameter E T = 500 GeV. 
We list the background process (BG), the range of hard scattering px {pt S ) over which they were 
evaluated, the number of events generated (N), number of events passing cuts (N cu t), total BG 
cross section (a to t in fb), and background cross section after cuts {a cut in fb). We compute the 
combined background to be a = 1.86 ± 0.36 fb, which yields a 95% CL upper limit of a = 2.44 fb, 
which we use for the computation of the LHC reach. We take m t = 170 GeV. 



BG 


FT 


N 


N cut 


^tot 


Pent 


QCD 


200-1000 


2.843 x 10 7 





8.9 x 10 7 fb 




QCD 


500-1000 


10 5 





1085 




QCD 


800-2000 


2 x 10 6 


17 


8.0 x 10 5 


(6.8 ± 1.6) x 10- 1 


QCD 


2000-3000 


2.5 x 10 5 


21 


103 


(8.7 ± 1.9) x lO" 3 


tt 


500-1000 


5.5 x 10 5 





7.2 x 10 5 




tt 


1000-2000 


1 x 10 5 


14 


0.48 


(6.7 ± 1.8) x 10" 3 


IL+jets 


10-800 


1.78 x 10 7 





3.1 x 10 7 




IL+jets 


300-800 


4 x 10 5 


9 


9.9 x 10 3 


(2.2 ±0.74) x lO^ 1 


IL+jets 


800-2000 


8 x 10 5 


1105 


0.89 


(1.2 ±0.04) x 10- 1 


Z+jets 


100-300 


5 x 10 4 





147 




Z+jets 


300-1000 


2.5 x 10 4 


6 


3 x 10 3 


(7.3 ±3.0) x lO" 1 


Z+jets 


1000-2000 


1 x 10 5 


1314 


7.0 


(9.1 ±0.25) x 10~ 2 
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TABLE III. Nine cases examined for S-multiplicity, and the possibility of reconstructing a 
Higgs mass via m(bb). We take Aq = 0, and fi < 0, except for case 4, for which /j, > 0. All mass 
quantities are in units of GeV. 



Case 


m 


mi/2 


tan (3 


mH e 


<n(6-tag)) 


(pr(6-tag)) 


(n (no-tag)) 


(pr(no-t 


1 


100 


300 


2 


87.7 


0.90 


127 


4.5 


124 


2 


100 


300 


10 


116.6 


0.90 


126 


4.5 


124 


3 


300 


300 


2 


88.4 


1.16 


121 


4.6 


121 


4 


300 


300 


2 


102.5 


1.18 


124 


4.7 


121 


5 


300 


300 


10 


116.6 


1.04 


123 


4.7 


121 


6 


600 


300 


2 


90.0 


1.65 


118 


5.1 


113 


7 


600 


300 


10 


117.2 


1.31 


124 


5.7 


119 


8 


1000 


275 


2 


92.3 


1.59 


126 


5.8 


123 


9 


1000 


350 


10 


120.5 


1.68 


128 


6.2 


122 
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FIGURES 



FIG. 1. Total cross section for strongly interacting sparticle pairs (solid), associated production 
of gluinos or squarks with charginos or neutralinos (dot-dashed) and chargino/neutralino produc- 
tion (dashes), for a) rriq = m g and b) rriq = 2m g , as a function of m g , for pp collisions at yfs = 14 
TeV. 

FIG. 2. Cross section in fb after cuts, as a function of cut parameter Ej,, for total background 
(x's), and for the six signal cases listed in Table 1. 

FIG. 3. a) Cross section in fb after cuts, as a function of Aq/tuq, for parameter choices shown. 
In b), we show variation in the cross section for Table 3 case 3 when tan/? is varied. 

FIG. 4. Background cross sections vs. Ej, obtained from various samples of events. The p^ s 
range of the hard scattering subprocesses is listed. 

FIG. 5. Signal cross sections for mo = (squares) and mo = 4m 1 / 2 (triangles), after cuts, 

using Ej, = 500 GeV, as a function of m g . We also show the 5a background cross section assuming 
an integrated luminosity of 10 fb -1 . We take Aq = 0, tan (3 = 2, and \i < 0. 

FIG. 6. In a), we show regions of SUGRA parameter space excluded by theory (bricks), ex- 
cluded by experiments (slashed), and regions explorable by LHC assuming integrated luminosity 
of 10 fb -1 , for Aq = 0, tan f3 = 2, and fi < 0. We also show regions explorable via slepton searches 
at LHC, and Higgs and chargino searches at LEP II. Much of the region below the spoiler mode 
(Z,2 — > Z\Z and Z2 — > Z\H{) contour is explorable via isolated trileptons at both Tevatron and 
LHC. In b), we show the same plane, but with various mass contours for m g and rriq (averaged 
over 1st generation). 

FIG. 7. Same as Figure 5, except for \i > 0. 

FIG. 8. Same as Figure 5, except for tan (3 = 10. 

FIG. 9. Fractional jet multiplicity after cuts (signal plus background), for six cases listed in 
Table 1, for different values of Ej,. The jet p? cut has been relaxed to PrQet) > 50 GeV. 

FIG. 10. Distribution in M est for different values of m g and cut parameter Ej,, along with 
background contribution (shaded) . M es t is formed after cuts by dividing the event into hemispheres 
via the transverse sphericity eigenvector, forming the invariant mass of the jets, and taking the 
maximal value. Events with only one jet in a hemisphere are rejected. We take ijiq = m^ 2 . 

FIG. 11. Same as Fig. 10, except we take mo = 4m 1 / 2 - 
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FIG. 12. Mass distributions for pairs of tagged b jets, (a) Case 1 from Table III. Higgs produc- 
tion is small, (b) Case 3 from Table III. B(Z2 — ► Z\Hi) = 94%, and peak is observed at mu e - (c) 
Case 8 from Table III. B(Z2 — > ZiHg) = 98%, but more combinatorial background. 
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